Двумерная адаптивная треугольная фильтрация в круглом окне с учетом районирования

Процедура позволяет разделять нестационарное двумерное поле на компоненты треугольным фильтром. В процессе обработки параметры окна фильтрации изменяются в соответствии с представлениями обработчика.

Алгоритм

Окно фильтрации имеет форму круга, диаметр которого задается извне и должен быть уже записан в некоторой ТОС. При корреляционном подходе к определению размеров окон эти данные могут быть получены процедурой расчета интервалов корреляции по ДАКФ (для перехода к диаметрам модуль интервалов корреляции необходимо домножить на 2 с помощью калькулятора свойств ТОС).

Конкретные весовые коэффициенты фильтра определяются размерами окна. Максимальное значение имеет коэффициент, расположенный в центре окна, минимальные - на краях. Соотношение максимального и минимального коэффициента - 10 к 1. Величина оставшихся коэффициентов изменяется линейно в зависимости от расстояния до центра окна. Полученные весовые коэффициенты нормируются так, чтобы их сумма была равна 1. Такая нормировка необходима для того, чтобы профильтрованное поле не отличалось по уровню от исходного.

При обработке больших массивов данных, контур которых включает территории с резко различными геологическими характеристиками, может возникнуть необходимость учитывать эти различия. При обработке с учетом районирования в окно фильтрации включаются только те точки, которые относятся к тому же району, что и центральная точка окна.

Такой подход заменяет следующую последовательность операций:

  • Разделение исходного поля на столько полей, сколько имеется районов, причем в каждом из новых полей оставлены значения исходного поля только в пределах одного района, а в остальных они заменены на пропуски;
  • Независимая обработка всех новых полей;
  • Синтез независимых результатов в единые результирующие поля.

Для каждой точки исходного поля производится свертка значений поля в скользящем адаптивном окне фильтрации с весовыми коэффициентами треугольного фильтра. Результат свертки записывается в Профильтрованное поле, разность между Исходным и Профильтрованным полем - в Остаточное.

Выходные значения могут быть неопределенными в двух случаях: когда значение исходного поля в центре окна не определено (пропуск) или когда в окне много пропусков (число пропущенных значений в три раза больше числа оставшихся).

Параметры

Для определения обрабатываемого поля необходимо выбрать нужное свойство из числа существующих в текущей ТОС (см. Рис. 1, параметр Свойство) на панели Исходное поле.

Рис. 1. Окно процедуры адаптивной треугольной фильтрации в круглом окне с учетом районирования

Для учета районирования обрабатываемой территории необходимо указать имя ТОС на панели Районирование, после чего выбрать нужное свойство из числа существующих в ней целых однобайтовых. Вообще говоря, учитывать районирование необязательно и для отказа от него достаточно снять отметку (флажок) об учете районирования, которая находится слева от поля ввода имени файла ТОС.

Диаметры окон должны быть определены до запуска процедуры и размещены либо в текущей ТОС, либо в другой ТОС, параметры сети которой совпадают с параметрами текущей. Панель Диаметры окон и позволяет указать источник информации.

Для определения Профильтрованного поля необходимо определить Имя файла ТОС и Свойство на панели Результат адаптивного осреднения в круглом окне с учетом районирования. По умолчанию используется имя файла текущей ТОС.

Остаточное поле определяется аналогично профильтрованному.

В зависимости от особенностей решаемой задачи, пользователь может отказаться от расчета одного из выходных полей: профильтрованного или остаточного - для этого достаточно снять отметку (флажок) о расчете соответствующей компоненты, которая находится слева от поля ввода имени файла ТОС.

В результате выполнения процедуры в зависимости от заданных пользователем имен выходных ТОС создаются новые ТОС, либо дописываются новые свойства к уже существующим ТОС, в которых размещаются профильтрованное и остаточное поля. Создаваемые ТОС наследуют размеры и параметры сети исходной (текущей) ТОС.